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Abstract 



The proliferation of different wireless access technologies, together with the growing number of 
multi-radio wireless devices suggest that the opportunistic utilization of multiple connections at the 
users can be an effective solution to the phenomenal growth of traffic demand in wireless networks. In 
this paper we consider the downlink of a wireless network with N Access Points (AP's) and M clients, 
where each client is connected to several out-of-band AP's, and requests delay-sensitive traffic (e.g., 
real-time video). We adopt the framework of Hou, Borkar, and Kumar, and study the maximum total 
timely throughput of the network, denoted by C T 3, which is the maximum average number of packets 
delivered successfully before their deadline. Solving this problem is challenging since even the number 
of different ways of assigning packets to the AP's is N M . We overcome the challenge by proposing 
a deterministic relaxation of the problem, which converts the problem to a network with deterministic 
delays in each link. We show that the additive gap between the capacity of the relaxed problem, denoted 



by Cdet, and C T 3 is bounded by 2^J N(Cdet + x)' which is asymptotically negligible compared to C^t, 
when the network is operating at high-throughput regime. In addition, our numerical results show that the 
actual gap between C T 3 and C^t is m most cases much less than the worst-case gap proven analytically. 
Moreover, using LP rounding methods we prove that the relaxed problem can be approximated within 
additive gap of N. We also extend the analytical results to the case of time-varying channel states and 
weighted total timely throughput. 
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I. Introduction 

Consumer demand for data services over wireless networks has increased dramatically in 
recent years, fueled both by the success of online video streaming and popularity of video- 
friendly mobile devices like smartphones and tablets. This confluence of trends is expected 
to continue and lead to several fold increase in traffic over wireless networks by 2015, the 
majority of which is expected to be video flU. As a result, one of the most pressing challenges 
in wireless networks is to find effective ways to provide high volume of top quality video traffic 
to smartphone users. 

With the evolution of wireless networks towards heterogeneous architectures, including wire- 
less relays and femtocells, and growing number of smart devices that can connect to several 
wireless technologies (e.g. 3G and WiFi), it is promising that the opportunistic utilization of 
heterogeneous networks (where available) can be one of the key solutions to help cope with the 
phenomenal growth of video demand over wireless networks. This motivates two fundamental 
questions: first, how much is the ultimate capacity gain from opportunistic utilization of network 
heterogeneity for delay-sensitive traffic? and second, what are the optimal policies that exploit 
network heterogeneity for delivery of delay- sensitive traffic? 

In this paper, we study these questions in the downlink of a heterogeneous wireless network 
with N Access Points (AP's) and M clients, where each client can request data from a subset 
of AP's. We assume that each AP is using a distinct frequency band. Therefore, the channels 
from AP's to the clients are orthogonal. Moreover, we assume that all AP's are connected to 



each other through a Backhaul Network (see Fig. la), with error free links, so that we can focus 
on the wireless segment of the problem. We model the wireless channels from every AP to each 
receiver as i.i.d. packet erasure channels. 

We focus on real-time video streaming applications, such as video-on-demand, video con- 
ferencing, and IPTV (TV distributed over an IP network), that require tight guarantees on the 
timely delivery of the packets. In particular, the packets for such applications have strict-per- 
packet deadline; and if a packet is not delivered successfully by its deadline, it will not be 
useful anymore. As a result, we focus on the notion of timely throughput, proposed in J3, 
which measures the long-term average number of "successful deliveries" (i.e., the packets that 
are delivered before the deadline) for each client as an analytical metric for evaluating both 
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throughput and QoS for delay-constrained flows. 

In this framework, the time is slotted and time-slots are grouped to form intervals of length 
r. For each interval every client has a packet to receive and the AP's have to decide on a 
scheduling policy to deliver the packets. If a packet is not delivered by the end of that interval, 
it gets dropped by the AP's. Timely throughput of a client is defined as the long-term average 
number of successful deliveries for that client. Furthermore, the total timely throughput (T ) 
is the summation of timely throughputs of all the M clients in the network. Our objective is 
then to find the maximum achievable T 3 , which we denote by C T 3 , over all possible scheduling 
policies. 

The challenge is that for each interval, even the number of different ways of assigning packets 
to the AP's is N M , which grows exponentially in the number of clients (M). Therefore, we face 
a quite complicated combinatorial optimization problem to find C T 3. To overcome this challenge, 
we propose a deterministic relaxation of the problem, which is based on converting the problem 



to a network with deterministic delays for each link. As we will show in Section III the relaxed 
problem can be viewed as a packing problem, in which each AP turns into a bin with certain 
capacity and each packet turns into an object which has different sizes at different bins. The 
relaxed problem is then to maximize the total number of objects that can be packed in the bins, 
denoted by Cd et . 

Our main contribution in this paper is two-fold. First, we prove that the gap between the 
solutions to the original problem and its relaxed version (i.e., the gap between C T 3 and C det is 
at most 2^J N(Cd S t + By noting that the number of AP's, N, is typically very small (in 
most cases between 2-4), the above result indicates that C det is asymptotically equal to C T 3 as 
C T 3 — > oo. Furthermore, our numerical results demonstrate that the gap between the proposed 
relaxed problem and the original problem is in most cases much smaller than the worst-case gap 
that we prove analytically. Therefore, instead of solving our main maximization problem we can 
solve its relxed version, and still get a value which is very close to the optimum. 

Second, we prove that the relaxed problem can be approximated in polynomial-time (with 
additive gap of N) using a simple LP rounding method. This approximation is appealing as 
N is usually limited and negligible compared to C^t- As a result, the solution to the relaxed 
problem provides a scheduling policy that provably achieves a T 3 that is within additive gap of 
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N + 2^N(C J3 -^) ofC j3 for C j3 > 7 -f. 

We also extend our problem formulation to the case of having time-varying channels and 
real-time traffic, where at the beginning of each interval clients have request for variable number 
of packets. We show that the aforementioned results hold in these two extensions, too. Moreover, 
when we allow for different flows to have different weights, and for bounded ratio of the weights, 
denoted by 0J max , we show that the gap between the weighted total timely throughput, C W _ T 3, 
and the solution to the relaxed problem, C w . det , is at most 2uj max \J N(C w . Aei + ^). 

Related Work: Although there are classical results 0, on scheduling clients over time- 
varying channels and characterizing the average delay of service, in recent years there has been 
increasing research on serving delay- sensitive traffic over wireless networks. This increase is 
due to the phenomenal increase in the volume of delay- sensitive traffic, such as video traffic. 
Dua et. al. [0 have considered weights and strict deadlines for packets; and if a packet is not 
delivered by its deadline, it causes a certain distortion equal to its weight. They have studied the 
problem of minimizing the total distortion and have characterized the optimal control. Agarwal 
and Puri [8J have considered a packet switched network where clients can get different types 
of service based on the amount they are willing to pay. Neely [9] has considered the problem 
of optimizing time averages in systems with i.i.d behavior over renewal frames. An algorithm 
which minimizes drift-plus -penalty ratio is developed in ||9). Shakottai and Srikant ifTOl have 
focused on minimizing the total number of expired packets, and have provided analytical results 
on scheduling. 

However, the most related work to this paper is the work of Hou et. al in [0 in 2009, in 
which they have proposed a framework for jointly addressing delay, delivery ratio, and channel 
reliability. For a network with one AP and N clients, the timely throughput region for the set of 
N clients has been fully characterized in [2] ; and the work has been extended to variable-bit-rate 
applications in [3], and time- varying channels and rate adaptation in Q. Although in ED- flU 
they provide tractable analytical results and low-complexity scheduling policies, the analyses are 
done for the case of having only one AP. This paper aims to extend the results to the case of 
having multiple AP's, where there is an additional challenge of how to split the packets among 
different AP's. 

The rest of the paper is organized as follows. Section|II]defines the notion of timely throughput, 
describes our network model, and explains the problem formulation. Section [III] describes and 
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formulates the relaxation of the problem, and states the main results. Section IV and Section [V] 



provide the proof of Theorem [T] and Theorem [2j respectively. Extensions to the main problem for- 
mulation are covered in Section |VTj Section [VrT| contains numerical analysis of our deterministic 



relaxation scheme. Finally, Section VIII concludes the paper. 



ii. Network Model and Problem Formulation 

In this section we describe our network model and precisely describe the notion of timely 
throughput introduced in 0. Finally, we formulate our problem. 

A. Network Model and Notion of Timely Throughput 

We consider the downlink of a network with M wireless clients, denoted by Rx 1; Rx 2 , . . ., 
Rxm, that have packet requests, and N Access Points APi, AP 2 , . . . , AP^r. These AP's have 
error- free links to the Backhaul Network (see Fig.l). In addition, we assume that time is slotted, 
and transmissions occur during time-slots. Furthermore, the time-slots are grouped into intervals 
of length r, where the first interval contains the first r time-slots, the second interval contains 
the second r time-slots, and so on. Moreover, each AP may make one packet transmission in 
each time-slot. 

Each AP is connected via unreliable wireless links to a subset (possibly all) of the wireless 
clients. These unreliable links are modeled as packet erasure channels that, for now, are assumed 
to be i.i.d over time, and have fixed success probabilities for each transmission. In addition, each 



channel is independent of other channels in the network. (In Section VI these assumptions will 
be relaxed to consider more general scenarios). The success probability of the channel between 
APj and RXj is denoted by py, which is the probability of successful delivery of the packet 
of RXj when transmitted by AP, during a time-slot. If there is no link between an AP and a 
client, we consider the success probability of the corresponding channel to be 0. Moreover, we 
assume that the channels do not have interference with each other. 

For now we assume that at the beginning of each interval each client has request for a new 
packet. Right before the start of an interval, each requested packet for that interval is assigned to 
one of the AP's to be transmitted to its corresponding client. Furthermore, during each time- slot 
of an interval, each AP picks one of the packets assigned to it to transmit. At the end of that 
time-slot the AP will know if the packet has been successfully delivered or not. If the packet 
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(a) (b) 



Fig. 1: Illustration of Our Network Model. Network configuration consisting of iV Access points 
(AP's), M wireless clients, packet erasure channels from AP's to the clients, and the Backhaul 
network is illustrated in (a). Our model for time, in which time is slotted and time-slots are 
grouped to form intervals of length r, is shown in (b). In this figure r = 4. 

is successfully delivered, the AP removes that packet from its buffer and does not attempt to 
transmit it any more. The packets that are not delivered by the end of the interval are dropped 
from the AP's. 

Definition 1: The decisions on how to assign the requested packets for an interval to the AP's 
before the start of that interval, and which packet to transmit on a time- slot by each AP are 
specified by a scheduling policy. A scheduling policy r\ makes the decisions causally based on 
the entire past history of events up to the point of decision-making. We denote the set of all 
possible scheduling policies by S. 

Definition 2: A static scheduling policy, denoted by 77 sta tic> is a scheduling policy in which 
each AP becomes responsible for serving packets of a fixed subset of clients for all intervals; 
and the packets of clients assigned to an AP are served according to a fixed order. In particular, 
a static scheduling policy 77 static is fully specified by a pair (II, T), in which II = [X 1 ,X 2 , . . . ,In], 
where Zj's partition the set {1,2,..., M}, indicating how the packet of clients are assigned to 
AP's. Furthermore, T specifies the ordering for the packets assigned to each AP. When ?7 sta tic is 
implemented, each AP is responsible for serving packet of the clients assigned to it by IT; and 
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each AP persistently transmits a packet until it is delivered successfully, before moving on to 
the packet of the client with the immediate lower rank in the ordering specified by T. 

Definition 3: A static scheduling policy is called greedy, and denoted by 77 g . s tatic> if the order 
of clients specified by T is according to the success probabilities of channels from AP to those 
clients, in decreasing order. 

Assume that a particular scheduling policy 77 is chosen. For any interval r (r G N), let 
N(r, 77) = [Ni(r, 77), A^(r, 77), . . . , Nm{t-, 77)] denote the vector of M binary elements whose j th 
element Nj(r,r)) is 1 if client RXj has successfully received a packet during the r th interval, 
and otherwise. When using scheduling policy 77, the timely throughput of RXj, denoted by 
Rj(i]), is defined as 

... V: . NJk.n) □ 

(1) 



RJri) 4 liminf ^k=i N i( k >v) ^ j = \ 2,...,M. 

r^oo r 



In simpler words, Rj{rf) is the long-term average number of successful deliveries for the j client. 
Further, we denote the vector of all i?j(r/)'s by R(r/), where we have R(r]) = [Ri(r]), i?2(^), • • • 
, Rm{v))- Therefore, the capacity region for timely throughput of M clients in the network is 
defined as follows: 

C 4 {^(77) : r] G S}. (2) 
Finally, the total timely throughput resulting from using 77, T 3 (r]), is defined as 

M 

T\v) = \\R(v)\\i = Y, R ^- ( 3 ) 
Therefore, T 3 (rj) denotes the timely throughput of the entire network. 

B. Main Problem 

Our objective is to find the maximum achievable total timely throughput, denoted by C T 3. 
More precisely, our optimization problem is 

Main Problem (MP): C j3 = sup 1 1 R(rj) 1 1 1 , (4) 



where ||i?(r/)||i = ^7=1 Rj(v)- Later in Section VI we will consider the problem of finding the 



maximum weighted total timely throughput J^jLi^jRjiv) an d i ts corresponding policy 77; but 
for now we focus on the problem in the case that u% = LO2 = ■ ■ ■ = % = 1. 

'More precisely, Rj{rf) = sup R s.t. limirifV^oo ^ fc=1 ^ 3< -* ; '' ? ^ > R with probability one. 
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As we state later in Lemma [T] in Section IV C T 3 can be achieved using a greedy static 



scheduling policy. Therefore, without loss of generality, the optimization in ([4]) can be limited to 
optimization over greedy static scheduling policies. So, the optimization in ([4]) is equivalent to 
finding the partition 11 such that the corresponding 77g. sta tic maximizes ||-R(% s tatic)||i- However, 
this is still quite challenging due to the fact that the number of possible greedy static scheduling 
policies to consider is N M , which grows exponentially with M. 

To overcome the aforementioned challenge, we propose a deterministic relaxation of the 
problem, which converts the problem to a network with deterministic delays in each link. We will 
show that the relaxed problem can be solved effectively, and its solution is asymptotically the 
same as the solution to our original problem. In the next section we will explain our deterministic 
relaxation approach, and then present the main result. 

in. Deterministic Relaxation and Statement of Main Results 

In this section we first explain the intuition behind our relaxation scheme and formulate the 
relaxed problem. Then, we state the main results. 

A. Deterministic Relaxation 

Consider the system model again: we assumed channel success probability Pij between APj 
and RXj, i = 1, 2, . . . , N, j = 1, 2, . . . , M. For now, suppose that APj wants to transmit a packet 
to client j, and that packet is the only packet assigned to APj for transmission. Furthermore, 
suppose that there is no deadline for transmitting the packet. Thus, APj persistently sends that 
packet to client j until the packet goes through. The number of time-slots expended for this packet 
to be delivered is a Geometric random variable where Pr(Gjj = k) = —pij) h ~ 1 , k = 
1,2,.... We know that ElGu] = — , and without any deadline, it takes — time-slots on average 

J Pij Pij 

for packet of RXj to be delivered when transmitted by APj. 

Therefore, a memory-less erasure channel with success probability can be viewed as a 
pipe with variable delay which takes a packet from APj and gives it to HXj according to that 
variable delay. The probability distribution of the delay is Geometric with parameter p^. 

To simplify the problem, we proposed to relax each channel into a bit pipe with deterministic 
delay equal to the inverse of its success probability. Therefore, for any packet of Poc,, when 
assigned to APj for transmission, we associate a fixed size of ^- to that packet. This means that 
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each packet assigned to an AP can be viewed as an object with a size, where the size varies 
from one AP to another; because ^-'s for different z's are not necessarily the same. On the 
other hand, we know that each AP has r time-slots during each interval to send the packets that 
are assigned to it. Therefore, we can view each AP during each interval as a bin of capacity 
r. Therefore, our new problem is a packing problem; i.e., we want to see over all different 
assignments of objects to bins what the maximum number of objects is that we can fit in those 
iV bins of capacity r. We denote this maximum possible number of packed objects by C& et . 
More precisely, if we define Xij as the — 1 variable which equals 1 if packet of client j is 
assigned to AP i5 and otherwise, then the relaxed problem can be formulated as following. 

N M 

Relaxed Problem (RP): Cdet = max Xjj 

i=l j=l 
M 

s.t. \"^i<T t = l,2,...,N (5) 

N 

X)zii<l J = 1,2,...,M 

i=i 

6 {0, 1}. 

B. Main Results 

We now present the main results of the paper via two Theorems. Theorem [T] gives bounds 
on the gap between the solution to the main problem Q and its relaxed version Q, and shows 
that one does not lose much if they solve the relaxed problem instead of the main problem. 
Furthermore, Theorem [2] proposes a performance guarantee for the polynomial-time algorithm 
which rounds the LP relaxation of the relaxed problem. The proofs of the two Theorems are 
provided in Section IV and Section |V} 



Theorem 1: Let C T 3 denote the value of the solution to our main problem in ([4]). Also, let 
Cdet denote the value of the solution to our relaxed problem in ([5]). We have 

C d et - 2^N(C det + j) < C v < C det + N. (6) 

Remark 1: The right part of the inequality in ([6]) suggests that C T 3 — Cdet can be no more 
than N. But the number of AP's N is limited and is usually around 2,3, or 4. Therefore, as 
Cdet — > oo ^=r- — > 0. Moreover, the left inequality in Theorem jlj suggests that C de t — C T 3 can 
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be no more than 2 J N(C det + =j). This additive gap becomes negligible compared to C det as 
Cdet — > oo. In addition, the inequalities in Theorem [T] imply that as C T 3 — >■ oo, Cdet — > oo, too. 
Therefore, ^ — > 1, as C T 3 — >■ oo. Hence, the bounds in Theorem |Tj suggest the asymptotic 
optimality of solving Cd e t instead of C T 3. 

Theorem [T] basically bounds the gap between C T 3 and Ca et . However, a remaining question 
is: if we run the system based on the greedy static scheduling policy which uses the assignment 
proposed by the solution to the relaxed problem, how much do we lose in terms of total timely 
throughput compared to C T 3? The following corollary which is proved in Appendix [E] addresses 
this question. 

Corollary 1: Assume C T 3 > Let n de t denote the assignment of clients to AP's suggested 
by the solution to the relaxed problem ([5]), and 77^ be the corresponding greedy static scheduling 
policy. Then, we have 



3iV 

Cj3 - N - 2^N(C J3 --£■)< H^deOlli < C T 3. (7) 

Remark 2: As we prove in Appendix [A] the upper bound given in the right inequality in 
Theorem [T] is tight. Furthermore, the lower bound given in the left inequality of Theorem 1 is 
tight in terms of order, i.e., there exists a network configuration and a positive constant k for 
which Cdet — C T 3 > ky/NCfet- 

The proof of tightness of the bounds given by Theorem [T] can be found in Appendix |A} 

Remark 3: The bounds in Theorem 1 are worst-case bounds, and via numerical analysis we 
observe that the gap between the original problem and its relaxed version is in most cases 
much smaller. Therefore, the solution to the relaxed problem tracks the solution to the main 
problem very well, even for a limited number of clients. To illustrate this, consider the network 



configuration in Figure 2a In this network there are two AP's with the same coverage radius |, 
and 10 clients which are uniformly and randomly located in the coverage area of the two AP's. 
The erasure probability of the channel between a client and an AP is proportional to the distance. 
Furthermore, the packets' deadline, r, is assumed to be 15. For 30 different realizations of this 



network, C T 3 and C det have been calculated, and plotted in Figure 2b (more detailed numerical 



results are provided in Section VII). 



The numerical results suggest that even for small-scale networks Cdet is usually very close 
to C T 3 . Therefore, rather than solving the main problem, we can solve the relaxed problem and 
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1 ' 1 1 1 ' 

5 10 15 20 25 30 

Enumeration of different realizations of the network 

W~ (b) 



Fig. 2: Numerical analysis for the gap between C T 3 and C det for the case of two AP's with 
coverage radius |, 10 wireless clients, and intervals of length r = 15. (a) illustrates the network 
configuration, where erasure probability of a channel is proportional to the distance between 
the AP and the corresponding receiver, (b) demonstrates the numerical results for the gap for 
30 different realizations of the network, where each realization is constructed from a random 
and uniform location of clients in the network. Each '+' indicates the value of C T 3 for each 
realization, while 'o' indicates the value of C det for the same realization. 



still get a value as the solution which is very close to the solution for the main problem. 

So far, we have shown by Theorem [T] that by considering the relaxed problem (RP) we do not 
lose much in terms of total timely throughput capacity. Nevertheless, in order for the relaxation to 
be useful there should be a way to solve the relaxed problem efficiently. The next Theorem, which 
is proved in Section |V} demonstrates that the relaxed problem can be approximated efficiently. 

Theorem 2: Suppose that x* tr is a basic optimal solution to the LP relaxation of We have 

N M 
i=l 0=1 

Remark 4: According to Theorem [2] if we find a basic optimal solution to LP relaxation of 
([5]), and then round down the solution to get integral values, the result will deviate from the 
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optimal solution by at most N. Note that N is typically very small (in most cases between 2-4); 
and therefore, this algorithm performs well in approximating the optimal solution of the relaxed 
problem ([5]). 

Remark 5: The relaxed problem in ([5]) is a special case of the well-known Maximum Gen- 
eralized Assignment Problem (GAP). There is a large body of literature on GAP; and its 
special cases capture many combinatorial optimization problems, and have several applications 
in computer science and operations research. GAP is APX-hard lfl4l . meaning that there is 
no polynomial-time algorithm with constant-gap performance guarantee for it. However, there 
are several approximation algorithms for GAP, including lfT4l . [fT51 . In particular, Chekuri and 
Khanna lfl4l . based on a modification of the work by Shmoys and Tardos lfT3ll . have come up 
with a 2-approximation algorithm for GAP; and Fleischer et. al. [15J have proposed an LP- 
based ^--approximation algorithm. The performance guarantees in the literature are concerned 
with multiplicative gap. However, our result in Theorem [2] suggests an additive gap performance 
guarantee of iV for the special case of GAP presented in ([5]). Since N (the number of access 
points) is typically very small, this provides a tighter approximation guarantee for our problem 
of interest. 

iv. Analysis of Approximation Gap (Proof of Theorem p]) 

In this section we prove Theorem [T] To this aim, we first state the following lemma which is 
proved in Appendix |Bj 

Lemma 1: C T 3 can be achieved using a greedy static scheduling policy. 
We now use Lemma [T] in order to prove the right part of the inequality in Theorem [T] 

A. Proof of C T 3 < C det + N: 

By Lemma [T] it is sufficient to prove that for any greedy static scheduling policy ^-static we 
have 

T 3 (?7g. s tatic) < Cdet + N. (9) 

Suppose an arbitrary greedy static scheduling policy 77 g _ stat i c with the corresponding partition 
fig-static = Pi,2^2, • • • ,%n] an d ordering r g . static is implemented. By ^ we know that 

M 

T 3 (r/g-static )=y~]Rj (%static ) • (10) 
3=1 
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On the other hand, by ([T]) we know that 

Rj (%stadc) = lim inf ELl Ni ( K %static) , jG{l,2,...,M}. (11) 

r— s>oo r 

Let Yi denote the random variable indicating the number of successful deliveries by APj during 
one interval, when 77 g - st atic is implemented, i = 1, 2, . . . , N. More precisely, 

Yi = Y N ^ ^- static )' * = 1> 2, . . . , iV. (12) 

Since a greedy static scheduling policy is implemented and channels are i.i.d over time, by LLN 
we have 

E[Y t ] = Vliminf U^&Ik^l = V^( %s ta ti c), 1 = 1,2,...,^. (13) 
Therefore, 

N N M 

Y E ^ = YY ^(%static) = ^(^g-static) = T 3 (r7g. static ). (14) 

i=l i=l jeJ 4 3=1 

Define ^ = |Z;|, and denote the enumeration of clients assigned to AP^ by {Ij(l),Xj(2), . . . ,Xi(qi)}, 
where the enumeration is according to the channel success probabilities of different clients in Zj. 
Let Gij be a geometric random variable with parameter pij, i — 1, 2, . . . , N, j — 1,2, ... , M. 
Then, it is easy to see that 

k 

Yi = max k s.t. y] Gg^) < r, i G {1, 2, . . . , iV}, k < q h (15) 

3=1 

since ?7 g - st atic persistently sends a packet until it is delivered, or the interval is over. Let us also 
define 

i 

h = max I s.t. Y 1 /PHi(j) < r > 1 < ft- ( 16 ) 

3=1 

Therefore, li is the maximum number of objects that fit into bin of capacity r when the channels 
are relaxed and clients in Zj are assigned to AP^. The following lemma (for which the proof is 
provided in Appendix [C]) relates l t to Y { . 

Lemma 2: Let r G N and G\, G2, ■ ■ ■ , G q be independent geometric random variables with 
parameters pi,p 2 , . . . ,p q respectively, such that 1 > pi > p 2 > . . . > p q > 0. Also define 

1 

Z = maxZ s.t. ^^l/pi<r (17) 

i=l 
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and 

i 

Y = maxi s.t. ^2Gj<t, i <E {1,2, ...,q}. (18) 

3=1 

Then, we have 

E[Y]<1 + 1. (19) 

Hence, by Lemma [2] we have 

N N 

J2 E i Y i\ + (20) 



Therefore, by (14) and (90) we have 

N 



T 3 (r/ g . sta tic) < J2 li + N - W 

i=l 

But, Yl!i=i h * s me value of the objective function in ^ for a feasible solution. Therefore, 

N 

J><C det . (22) 



i=l 



Putting (91 ) and (92) together we get 

T 3 (?7g. s tatic) < Cdet + N. (23) 

Hence the proof of the right inequality in Theorem [T] is complete. 



B. Proof of C det - 2^N{C det + f ) < C T s 

Consider the assignment proposed by the solution to the relaxed problem in ([5]), where 
the clients that are not assigned to any AP for transmission are assigned to AP's arbitrarily. 
Let Ilg^tatic = [Zf^Z^ 61 , . . . ,Z^ £ ] denote the resulting partition, and also let 7?g_ e s tatic denote the 
corresponding greedy static scheduling policy. Therefore, we have 

T 3 (^ atlc ) < C T3 . 



So, it is sufficient to prove that C det -2y / N(C det + f ) < T 3 (^ ta tic)- Let Y i det denote the random 
variable indicating the number of successful deliveries by APj during one interval, when ??g_ e s tatic 
is implemented, i — 1, 2, . . . , N. With the same argument as in part A we have 

TV 

T 3 «l tic ) = E^K detl 

i=i 
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Therefore, it is sufficient to prove that 

I N 

C det -2J N(C det + -) < J>K det ]. 

i=i 

Define g» = |Xf et |; and denote the enumeration of clients assigned to APj by {Xj(l),Xj(2), . . . 
,Xj(gj)}, where the enumeration is according to the channel success probabilities of differ- 
ent clients in Xj. Further, let Gy be a geometric random variable with parameter pij, i = 

1, 2, . . . , N, j — 1, 2, . . . , M. Then, it is easy to see that 

k 

y 4 det = max k s.t. ^ G iXi(jY « < r, i G {1, 2, . . . , N}, k < q h (24) 

3=1 

since ^™. static persistently sends a packet until it is delivered, or the interval is over. Also define 

t 

/ det = max I s.t. < r, Z < q-j. 

Therefore, Z det is the maximum number of objects that fit into a bin of capacity r when the 
channels are relaxed and clients in X det are assigned to APj. The following lemma (which is 
proved in Appendix [p]) relates Z det to Yj det . 

Lemma 3: Let r G N and Gi, G 2 , . . . , G g be independent geometric random variables with 
parameters p 1 ,p 2 , . . . ,p g respectively, such that 1 > p t > p 2 > . . . > p q > 0. Also define 

i 

I = max/ s.t. < r 

i=l 

and 

i 

Y = maxi s.t. Gj < r, i £ {1,2,..., q}. 

Then, we have 

1-2^1+^ <E[Y}. (25) 

By Lemma [3] we have 

AT N j N 

E /det - 2 E \A det + 4 < E^K det ]- < 26 ) 

i=l i=l i=l 

According to Cauchy-Schwarz ineqality we have 



i=l 



\ 



N N 
i=l 
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On the other hand, note that 

N 



C det = ^/, det . (28) 



i=i 



Putting (26), (27), and (28) together we get 



N 

C det -2yiV(C det +-)<C T 3. (29) 
Hence, the left inequality of Theorem [T] is proved and the proof of Theorem [T] is complete. 



v. Efficient Solution To The Relaxed Problem 
(Proof of Theorem [2]) 

In this section we show that the approximation proposed by rounding down a basic optimal 
solution to LP-relaxation of RP in ([5]) is at most N away from C det . Note that RP is a mixed 
integer linear program. Linear relaxation of such a problem replaces the integrality constraints 
with non-negativity constraints; i.e., it replaces the constraint Xij G {0, 1} with x# > for 
i — 1, 2, . . . , N, j — 1,2, ... , M. Let LR-RP denote the linear relaxation of RP. Any feasible 
solution to LR-RP can be denoted by an N-by-M matrix x = [xjj]jv*M sucn that each element 
of x, x^, is the fraction of the j th object which is assigned to the i th bin; if = 1, then the j th 
object is fully fit into the i th bin; but if is a fractional value, then the j th object is partially 
assigned to the i th bin. 

Note that finding a basic optimal solution to a linear program is straightforward, and is 
discussed in |fT8l . So, let x* tr = [x*j T ]N*M denote a basic optimal solution to LR-RP with 
objective value V*; i.e., V* = J2f=i ^2j=i x *if '■ Define K as the number of objects that are not 
assigned to any of the bins (not even partially), by the solution proposed by x* tr . Therefore, 
V* < M - K; and by noting that C det < V* (due to optimality of x* tr ), we have 

C det <M-K. (30) 

Moreover, suppose that there are Z objects in the solution proposed by x* tr each of which is 
associated with at least one fractional value x*j r . Therefore, 

N M 

EEL<fJ =M-K-Z. (31) 
i=i j=i 
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By putting ( [30] ) and ( |3"T| ) together we get 

N M 

^et-EE^J < Z ( 32 ) 

t=l j=l 

Now, it is sufficient to show that Z < N in order to complete the proof of Theorem [2j To do 
so, we convert RP to an instance of Min-GAP, which is formulated as following, and then use 
the result in |fl6l iflTI 

N M 



Min-GAP: min 

i=l j=l 



M 



s.t. SijXij < b i — 1, 2, . . . , N (33) 

N 

J2 X H = 1 J = 1,2,...,M 

i=i 

ar y G {0, 1} i = 1, 2, . . . , N and j = 1, 2, . . . , M. 

The conversion is done as following; and it is similar to that of Section 3.2 in [fT4ll : 

• Set Sij = jr if p^- > 0, and Sjj = T + 1 if p^- = 0. 
. Set b = t. 

. Set dj = 1 for i = 1, 2, . . . , N and j = 1, 2, . . . , M. 

• Add an additional bin of capacity r to ensure the existence of a feasible solution. For this 
bin, the (N + l) st bin, set s N+lj = and c N+lj = M for j = 1, 2, . . . , M. 

Therefore, the converted problem is as follows. 

Converted Problem (CP): min 2_. x ?i + ^"^(jv+iJj 

»=i 3=i i=i 

M 

s.t. ^SijXij <r i = 1,2, N + 1 (34) 

7 = 1 



5^a;y = l J = 1,2,...,M 



i=l 



Xij G {0, 1} z = 1, 2, . . . , JV + 1 and j = 1, 2, . . . , M. 



Let LR-CP denote the linear relaxation of the converted problem in (34). The following lemma, 
which is proved later in this section, connects basic optimal solutions of LR-RP and LR-CP. 
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Lemma 4: Suppose that x* tr = [x*j r ]7v*M is a basic solution to LR-RP. Further, consider 

x * = [ x *j}(N+i)*M, where 



1) Vi e {1, 2, . . . , iV} and Vj G {1, 2, ... , M} x% = x 



*tr 
ij ■ 



2) Vje{l,2,...,M} x* N+lj = l-J2lt 



N „*tr 



X; 



Then, x* is a basic solution to LR-CR 

Consider the solution x* to LR-CP with the same construction as described in Lemma |4j By 
Lemma |4j x* is a basic solution to LR-CP. Therefore, by considering the fact that there are also 
exactly Z partially assigned objects in the assignment proposed by x*, and using Theorem 1 of 
Ifl6t ifTTl we get 

Z < Number of bins that are full. (35) 



Since the (N + l) st bin is never used up to its capacity, by (35) we have 



Z < N. (36) 



Hence, by ( [32] ) and ( [36] ) we get 

N M 
i=l j=l 

Therefore, the proof of Theorem [2] is complete. We now prove Lemma [4] 

Proof: It is easy to check that x* is a feasible solution to LR-CP. Consider the feasible 
region for a^/s, 1 < i < N,l < j < M, in LR-RP and LR-CP; and denote the corresponding 
regions by IZi and 1Z 2 , respectively. Note that each point inside TZi corresponds to a feasible 
solution for LR-RP; and each point inside IZ2 corresponds to a feasible solution for the first N 
bins of LR-CP. It is then straightforward to check that TZi = 1Z 2 . In addition, any basic solution 
to LR-RP is an extreme point of the polytope defined by IZi. Now suppose that x* is not a basic 
solution to LR-CP. Then, there exists a constant < a < 1, and two non-identical matrices 
w = [?%-](7v + i)*Af,y = [yij)(N+i)*M such that 

x^ j = oew ij + (l-a)y ij , i — 1,2, . . . , N + 1 j = l,2,...,M. (38) 

Therefore, since x*f = x*j for 1 < i < N, 1 < j < M, x* tr can be written as a convex 
combination of two distinct points inside TZ 1 , which is contrary to the fact that x* tr is the 
extreme point of the polytope TZ\. Hence, the assumption that x* is not a basic solution to 
LR-CP is incorrect, and the proof is complete. ■ 
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Corollary 2: Suppose we choose a basic optimal solution to the LP relaxation of ([5]), denoted 
by x* tr , and round down the solution to get integral values. Let n^* denote the assignment 
suggested by the resulting integral values; and let rfl* denote the corresponding greedy static 
scheduling policy.. For C T 3 > we have 

C T3 -2N- 2^N(C T 3-™) < \\R(vZ)\\i < C v . (39) 

Proof: Let = 2~2f=i 2~2jLi l x ij T \ denote the resulting objective value of the rounded 
down basic optimal solution to LR-RP. According to Theorem [I] and Theorem [I| > C^t — 
N > C T 3 — 2N. Therefore, by using the similar argument as the one in Corollary [T] the proof 
will be complete. ■ 

vi. Extensions 

In this section we investigate two important extensions to our main problem formulation. First, 
we extend our model to allow for time-varying channels, where channels' statistics change over 
time. In addition, we consider real-time traffic pattern in which for each interval a client has 
request for a variable number of packets (possibly zero). Second, we investigate the problem of 
maximizing weighted total timely throughput. We present similar results as the main result in 
Theorem Q] for both extensions. 

A. Time-Varying Channels and Real-Time Traffic 

So far, we have considered a fixed packet generation for each interval, meaning that we have 
assumed that at the beginning of each interval each client has request for exactly one packet. 
This assumption can be modified by considering a time-varying packet generation pattern, in 
which for every interval, each client might have request for no packets, or it might have request 
for multiple packets. In addition, the number of packets requested by clients for one interval 
might depend on the number of packets requested for other intervals. 

Furthermore, so far we have assumed that channel reliabilities are static, and channel success 
probabilities do not change from time to time. But, this model can be generalized to include 
time- varying channels whose statistical behaviors are not necessarily independent of one another. 

We capture the above two generalizations by considering an irreducible Finite-State Markov 
Chain (FSMC), in which each state jointly specifies the subset of clients that have packet request, 
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and the number of their requested packets, as well as the channel states for different channels 
during an interval. When a new interval begins, the Markov Chain might change its state, and in 
this case, packets for a new subset of clients are requested, and the channel reliabilities change. 
By considering this FSMC we allow for variable-bit-rate traffic, as well as variations in channel 
reliabilities over time. 

Let us denote the set of all possible states of the FSMC by C. Each state A G C specifies a pair 
(5(A), P(A)), where 5(A) = [#i(A), B 2 (X), B M {\)\, such that B^X) is the number of the 
packets requested by client j, and P(A) is an N * M matrix that contains channel success 
probabilities. It is assumed that channel success probabilities remain the same during each 
interval, and are known to the AP's at the beginning of each interval. 

Our objective is again to find C T 3 . We use a similar argument as the one in [4] for extensions 
to time-varying channels and variable-bit-rate traffic. In particular, we decompose the set of 
intervals into different subsets, where each subset contains the intervals that are in the same 
state of the FSMC. 

Consider all the intervals for which the system is at state A. For those intervals we convert our 
problem to an instance of the problem described in Section |n} More particularly, for the system 
described by state A, we ignore all the clients that do not have packet request. Furthermore, for 
an Y j £ {1,2, . . . , M} where Bj(X) > 1 we consider £>j(A) — 1 virtual clients, such that the 
channel between APj and each of those virtual clients would have success probability Pij(X). 
This means that these virtual clients are exact copies of RXj. Consequently, for the intervals for 
which the system is at state A the problem becomes the same as the one described in Section [ITJ 
With the same argument as in proof of Theorem [IJ there exists a fixed assignment n(A), which 
if used together with its corresponding optimal ordering for such intervals, achieves the optimal 
T 3 for those intervals. We denote this optimal T 3 by C T 3(A). In addition, let Cd et (A) denote the 
solution to the relaxed problem when the system is at state A. For any state A G C, with the 
same argument as in the proof of Theorem [T] we have 




(40) 
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Now, let 7i\ denote the steady state probability of state A. Therefore, 

C T3 = 5> A C T3 (A) 

Aec 

Cdet = y^7r A C de t(A). 

AeC 

Hence, 



Cdet " 2 JiV(C det (A) + f ) < < ^det + N. (41) 



4 

Aec 



On the other hand, by using Cauchy-Schwarz inequality we have 



^vr A JiV(C det (A) + ^; 
Aec 



Aec V AeC 



N(C det + ^). (42) 



Putting (41 ) and (42) together we get 



N 

C det - 2^ N(C det + -) < C J3 < C det + N, (43) 

which is the same as the result in Theorem 1. 

Theorem 3: For the network model described in Section |n] consider the extension to time- 
varying channels and real-time traffic, modeled by the FSMC described in Section [VI| where 
each state of FSMC captures both the success probability of channels and the number of packets 
for each client during an interval. We have 



iV 

C det - 2yj N(C det + j)< C j3 < C det + N. (44) 

B. Weighted Total Timely Throughput 

In Section [TT] we considered the same importance for all the flows in the network; and our 
objective was to maximize T 3 . However, it might be the case that in a network some of the 
flows are more important than the others. Therefore, those flows should somehow be prioritized 
accordingly. 

In this section the formulation remains the same as the one described in Section |ll} except the 
objective function, which rather than maximizing T 3 , maximizes a weighted average of timely 
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throughputs. In particular, weighted total timely throughput of the scheduling policy rj, w-T 3 (r)), 
is defined as 

M 

w-J 3 (rf) = ^2ujRj(r]), (45) 
i=i 

where tu/s are arbitrary weights greater than 1 (j = 1,2,..., M). 
Our objective is to find 

C W _ T 3 = sup w-T 3 (r]). (46) 

For this extension of the problem we again propose the channel relaxation which results in 
a new integer program. This integer program is again a GAP. The formulation of the relaxed 
problem is as follows: 

N M 

C w -det = max ^ y ^ -Eij^j 

i=l j=l 
M 

s.t. y"^l< T i = l,2,...,JV (47) 

N 

^2xij<l j = l,2,...,M 

i=l 

x fi e{0,l} i = l,2,...,iV j = l,2,...,M. 
The following theorem, which is proved in Appendix |FJ states that the value of the solution 



to (46) is asymptotically the same as the value of the solution to (47) for C w _ T i — >■ oo (or 



equivalently C w .d et — > oo). 



Theorem 4: Let C W _ T 3 denote the value of the solution to (46). Further, let C w . det denote the 



value of the solution to (47). Then 



N 

C«i-det 2uj max 

where u max = maxjwi, u 2 , ■ ■ ■ , com}- 

vii. Numerical Analysis 

In this section we provide numerical analyses for our deterministic relaxation scheme. So, 
we consider a wireless network with two Access Points , and several wireless clients that are 
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10 





(c) (d) 

Fig. 3: Numerical Results, (a) illustrates the network configuration in which there are two AP's 
with coverage radius | each, M randomly and uniformly located clients in the coverage area of 
the AP's, and channel erasure probabilities proportional to the distances, (b) compares C T s with 
the T 3 which is resulted from the assignment proposed by ^ for the case of M = 10, r = 15 and 
10 different realizations of the network. '+' in (b) indicates the value of C T 3 for each realization, 
while 'o' indicates the value of T 3 (r] det ). (c) compares C^ t (denoted by '+') with the objective 
value of the rounded basic optimal solution (denoted by 'o') for the case of M = 20, r = 30 
and 20 different realizations of the network. Finally, (d) compares C T 3 (denoted by '+') with the 
T 3 resulted from the assignment proposed by rounding the basic optimal solution (denoted by 
'o') for the case of M = 10, r = 15 and 10 different realizations of the network. 
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uniformly and randomly located in the network (see Figure |2aJ). The channel from every AP 
to every wireless client is an erasure channel with erasure probability which is proportional to 
the distance between the AP and the client. The distances in the network are normalized, and 
we assume that the AP's have the same coverage radius R = |. Therefore, the channel erasure 
probability is 1 for the channel between an AP and a client which is located at the distance 
R = | from it. Furthermore, the distance between the two AP's is |. 

Figure |3b|corresponds to the case where M — 10 and r = 15. In each realization 10 clients are 
randomly located in the network. For each realization C T 3 is calculated. Then, the corresponding 
relaxed problem is solved, and the network is run for 10000 intervals under the assignments 



proposed by its deterministic relaxation solution. Fig 3b shows the comparision between the two 
for 30 different realizations of the network. 



Figure 3c demonstrates how our proposed LP-rounding algorithm performs compared to Cdet. 
We consider M = 20 and r = 30. 30 different realizations of network are considered. For each 
realization Cd e t is found. Then, we find the number of packets that are completely assigned to 
a bin in a basic optimal solution to the LP relaxation of ([5]). The result confirms the fact that 
our proposed algorithm performs well. The performance improves as we have larger number of 
clients in the network. 

Figure 3d shows how far our T 3 will be if we use our rounding approximation as the 



assignment strategy for the packets, and run the network for 10000 intervals according to that 
assignment. In this case we have considered 10 clients and have looked at 10 instances of the 
network, r is assumed to be 15 in this case. 



viii. Conclusion 

In this work we investigated the improvement by utilizing network heterogeneity in order to 
enhance the timely throughput of a wireless network. In particular, we studied the problem of 
maximizing total timely throughput of the downlink of a wireless network with N Access points 
and M clients, where each client might have access to several Access points. This problem is 
challenging to attack directly. However, we proposed a deterministic relaxation of the problem 
which is based on converting the problem to a network with deterministic delay for each link. 

First, we showed that the value of the solution to the relaxed problem, Cdet, is very close 
to the value of the solution to the original problem, C T 3. In fact, as C T 3 — >■ oo, ^ — >■ 1. 
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Furthermore, the numerical results indicate that for networks with limited number of clients, 
the gap between C T i and Cd et is very small. Second, we proposed a simple polynomial-time 
algorithm with additive performance guarantee of N for approximating the relaxed problem. This 
approximation performs well as the number of Access points is for most cases between 2-4. We 
also extended the formulation to allow time-varying channels and real-time traffic, and showed 
the same results for these extensions. Furthermore, for the problem of maximizing weighted 
total timely throughput we came up with similar results on the bounds on the gap between the 
original problem and its relaxed version. The future works will investigate more general network 
models, and application of our proposed deterministic relaxation to other frameworks . 
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Appendix A 

Proof of the tightness of the bounds in Theorem 1 

We first prove that the upper and lower bounds given in ([6]) are tight. More specifically, we 
show that there exist N, M, and some channel success probabilities for which C T 3 gets arbitrarily 
close to Cdet + N. In addition, there exist N, M, and some channel success probabilities for which 
0(|C T 3 — Cdet|) = 0(y/NCfet). Then, we prove Theorem 2. 

A. Proof of the tightness of the upper bound 

We show that for any given iV and < e < iV there exist M, r, and channel success 
probabilities such that C T 3 — C det = N — e. We set M = Nt, and we choose C det such that C7 det < 
M — N and y 6N. Further, for the channel between APj and RXj we set the channel success 
probability = Cia ^~ e , where % = 1, 2, . . . , iV and j = 1, 2, . . . , M. Therefore, according to 
symmetry, both the optimal assignment which results in C T 3 and the optimal assignment for the 
relaxed problem which results in C det assign r packets to each AP. Furthermore, without loss of 
generality we can assume that for APj packets of clients j = 1 + (i — l)r, . . . ,ir are assigned 
to APj. It is easy to check that the following inequalities hold for any APj, i — 1, 2, . . . , N: 

l+(i-l)r+CWJV 1 „ Ar l+(i-l)r+C dct /iV+l 

j=l+(?.-l)r J j=1+(i-1)t 

Therefore, the maximum number of packets that can be packed in the relaxed problem is C de t. 
Now, we calculate the expected number of packet deliveries: For any APj the expected number 
of successful deliveries during one interval is T~{ CiA ~^~ e ) = c ^ + ^~ e , Therefore, we have C T 3 = 

N( c ^ N N - e ) = C7 det + N-e. Hence, C j3 - C det = N - e. 
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B. Proof of the tightness of the order of the lower bound 

We show that there exists a wireless network realization for which O ( | C T 3 — Cdet \) — 0( V iVCdet) • 
More specifically, for a given N we show that there exist a positive constant k along with M, r, 
such that C det - C T 3 > ky/NC det . We choose C det such that ^GN, and we set M = C det . 
In addition, we set the channel success probability p i3 ■ = p = j& < 1 for some r e N, where 
i = 1, 2, . . . , TV and j = 1, 2, . . . , M. 

Therefore, both the optimal assignment which results in C T 3 and the optimal assignment for 
the relaxed problem which results in Cdet assign packets to each AP. It is easy to check that 
our chosen C det is actually the solution to the relaxed problem. Now, let Y denote the number 
of successful deliveries for one of the AP's. Thus, C T 3 = NE[Y]. Also, let I denote the number 
of packets that can be packed in a bin corresponding to a certain AP. Therefore, I — and 
p—\- We only need to show that there exists a constant k such that 



I - E[Y] > k\Tl. 



(50) 



Noting that / 



pr we have 



/ - E[Y] 



pr-Ei( 7 V(i-pP+' E C^^-pY'*] 

j=l j=l+l 

^ r E ( T liV" 1(1 " p)T " i_z ^ ^Va-pr - ' 

n± (;:>-'(i-#-'-E(( T ;V(;:!)«i-pn 

j=l+l KJ / 3=1+1 \ J ' \J / 



/[E JV^a-pr* 1 - E ( r ? V^-n 

j=z+i ^ ' 3=1+1 \ j y 




Now note that ( r l x ) 



^(I) = (l-p)(I). Therefore, 



i!(r-l-0! 
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By Theorem 2.6 of [19| we know that for positive integers m,n,q, with m > q > 1 and n > 1 

mn\ 1 _J_/ J__ i i_n _i m mn+\ 

^ m a m — a > 



J y/2% ' (m - q)( m -<i) n +hqi n +l' 



Substituting n by 1, m by r, and q by I we get: 



i 



T\ 1 T T 2 1/1 1 In 

> — - — e 12 t ! 



I J V2tt (r - Z) (t_z)+ ^ + 



2 



1 r^+f 



— f 1 - 1 - — ) 

_gl2\r I T-i ' 



2tt (r(l -p)) T - J+ 3(pr) ,+ 3 



However, - — T S > — 2. Therefore, 

' T I T — l 



All l _i 



Hence, we get 



I — E[Y] = jf T j 1 W-p) T -'=r^y +1 (l-pr 



> r ^ 1 1 e" V +1 (l - p) t ~' +1 

T vWVrp(l-p)^(l-P)^ P 1 PJ 



Thus, by setting = e e w ^-2 the proof will be complete. 

Appendix B 
Proof of LemmaQ] 

Given that clients are assigned to AP's according to some partition, the problem of finding 
the optimal scheduling policy is decomposed into finding the optimal scheduling policy for each 
of the AP's; since AP's are operating at different frequency bands, and have different packets to 
transmit. In [11] it has been shown that the timely throughput region for the clients assigned to 
an AP is a polymatroid. Thus, maximum of the summation of timely throughputs for the clients 
assigned to an AP is achieved at one of the extreme points. Furthermore, it has been shown 
in lfT2l that each extreme point of this polymatroid is achieved by a static scheduling policy. 
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Therefore, there is a static scheduling policy which achieves optimality. Note that when a static 
scheduling policy is implemented since the channel success probabilities are i.i.d. over intervals, 
by LLN the total timely throughput is equal to the expected number of packet deliveries during 
one interval. Therefore, the optimal scheduling policy is the static scheduling policy which results 
in the maximum expected number of delivered packets for one interval. 

We prove that for a given assignment n = [X 1 ,X 2 , . . . ,X N \ the optimal ordering of the packet 
of clients assigned to APj is according to the order of Pij, j E Xj. To prove this, it is sufficient 
to prove that for any given order if we swap the order of two adjacent clients in such a way 
that the client with the higher corresponding p^ is prioritized higher, then the expected number 
of deliveries will be no less than before swapping. 

Lemma 5: Let r £ N and G 1 ,G 2 , ■ ■ ■ ,G q be independent geometric random variables with 
parameters p%,P2, ■ ■ ■ ,Pq, respectively. Suppose that pa < pd+i for some d E {1, 2, . . . , q — 1}. 
In addition, let G[, G' 2 , . . . , G' q be independent geometric random variables (and independent of 
Gi's) with parameters p h p 2 , ■ ■ Pd-i,Pd+i,Pd,Pd+2, ■ ■ ■ ,p q , respectively. Then, 

i=l j=l i=l j=l 

i=l j=l j=l i=d+l j=l 

= EPr(EG;<r) + Pr(G, + EG;<r)+ E Pr (E^<-) 

i=l j=l j=l i=d+l j=l 

, > d—1 i d—1 q i 

< E Pr (E g 'j ^ r ) + Pr (^ + E G ^ r )+ E Pr (E G 'j ^ r ) 

i=l j=l j=l i=d+l j=l 

= E Pr (E G ^ r )' 

i=l j=l 

where (a) follows from the fact that success probability of Gd, which is pa, is less than success 
probability of G' d , which is Pd+i- ■ 
According to Lemma [5j the ordering which results in the maximum expected number of 
delivered packets is the one according to clients' channel success probabilities. So, when a 
specific partition is chosen, the optimal scheduling policy is the corresponding greedy static 
scheduling policy. 



Proof: We have 

i=l j=l 
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There are N different possible assignments for an interval, each when used together with 
its corresponding optimal ordering results in a certain expected number of delivered packets. 
Consequently, the optimal solution to the problem of which assignment (partition) to choose for 
each interval is simply to choose the assignment which results in the highest expected number 
of delivered packets, and to repeat it for all intervals. 

Appendix C 
Proof of Lemma[2] 

For I = 0, we have Pi < -, for % = 1, 2, . . . , q. Therefore, E[Y] in this case is less than that 
of the case in which pi — p 2 — . . . — p q — -. On the other hand, for pi = p 2 = . . . = p q = - 
E[Y] < t * - = 1. Hence, the statement is true for I = 0. Now, suppose that I > 0. We know 

that I = max I s.t. Y^i=i 1/Pi — T - Therefore, we have 

i n i+i n 

E 1 ^-<E 1 - (53) 
ttpi tip* 

We will show that E[Y] can be at most I + 1. Without loss of generality we can omit pj's that 
are equal to zero; because by omitting them neither of E[Y] nor / change, and E[Y] — I would 
remain the same. So, we suppose that 1 > p 1 > p 2 > . . . > p q > 0. It is sufficient to prove the 
lemma for the case of q = r; because if we have less than r geometric random variables, E[Y] 
will be less. On the other hand, we do not need to consider the case q > r; since for i > r, 
P r (X]j=i Gj — T ) = 0- Therefore, we suppose that q — r. 

Let Xi = Yl\=i Gi, where Gi = Geom(pi). By this notation we have: 

i-i 

Pr (X t > t) = Pr ( Y = i) ■ (54) 

i=0 

Now we write down the expression for E[Y): 

l-l T T 

E[Y] = J2 iFl ( Y = *) + E Pr ( X/ = 0(E iPr ( y = j \ Xl = *))■ (55) 

8=0 t=l 1=1 

Since 1 > pi > pi +1 > . . . > p T > 0, E[Y] is less than the case where p t = pi +1 — . . . — p T , 
although I remains the same. So it is sufficient to prove Theorem 1 for the case where p\ = 
pi + i = . . . = p T . For t < r if we set pi = pi + i — . . . — p T we have 

T 

^Pr(F = i\X t = t)= E[Y\Xi = t] = l + (r- t) Pl . (56) 

i=l 
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Therefore, by using ( [55] ) and ( [56] ) we have 

l-l T 



E[Y] = * Pr(^ = i) + +Pi( T ~ *)) Pr ( X < = *) 

i=0 t=i 
2—1 oo oo 

= ^zPr(y = i) + (l + Pl r)(l - Pr(X ; > r)) - Pl [£\Pr(X, = t) - ^ * Pr PG = *)] 

i=0 t=Z t=r+l 

2—1 2—1 2 . oo 

= ^iPr(y = i) + (Z+p I r)-(Z + p,r)^Pr(y = i)-p I ^-+p 1 ^ tPr(X / = t) 



i=0 i=0 i=l ^ ^ t=r+l 

2—1 2 oo 

= ^(i - I - pir) Pi(Y = i) + (I + Pl (r - -)) + Pi J2 1 Pr ^ = *) 

i=0 i=l P % t=r+l 

2—1 oo 

< ^(i-Z-p z r)Pr(K = i) + l + l+ Pl tPr{Xi = t), (57) 

8=0 t=T+l 

where the last inequality (a) follows from r < ~. an ^ the assumption that pi + i = pi. Now, 
we only need to rewrite p\ J2tL T +i t P r (^z = t) in terms of Y. For t > r we have 



2-1 

Pr(X, = t) = J2 Pl ( X i = A Y = i) Pr (^ = *)■ (58) 

i=0 

Therefore, 

oo oo 2— 1 

^tPr(X=t) = t(J]Pr(X, = t|y = z)Pr(y = i)) 

t=T+l t=T+l i=0 

2—1 oo 

= J^Pr(y = i)( tPi{Xi = t\Y = i)). (59) 

i=0 i=r+l 

But due to memory less property of geometric distribution, we know that 

oo oo 2 

^(t-r)Pr(X = t|y = z)= ^(t- r )Pr(^ Gj—t-r) 

t=T+l t=T+l j=i+l 

oo 2 2 1 

= E tPr (E G J = t ) = E Vi<Z-l. (60) 

Pi 



t=l j'=i+l i=i+l 



Therefore, 



Hence, 



tPr(X l = t\Y = i)=r+ -■ ( 61 ) 



t=T+l j=i+ 
oo 2—1 2 



^ tPr(X = t) = $^Pr(y = i){ Pl r+ -)■ ( 62 ) 



t=r+l i=0 J=i+1 
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Substituting ([62]) into ([57]) we get 



i-i i-i i 

E[Y] < ^(i-Z- Pl r)Pr(y = z) + / + 1 + ^ Pr(F = i)(prr + E — ) 

i=0 j=0 j=i+l ^ 

2-1 2 

= / + l + ^Pr(F = i){i-l-p l r + p l r + E -) <i + l, (63) 

i=0 3=i+l J ' 

where the last inequality follows from the fact that Vj G {i + 1, . . . , 1} Pi <Pj. 

Appendix D 
Proof of Lemma[3] 

We will show that E[Y] is more than / — 2yjl + |. It is sufficient to prove Lemma j^j for the 
case of q = /; because for q > I, E[Y] would only increase. On the other hand, q cannot be less 
than I according to the assumption I = max/ s.t. Yl\=i Vp» < T - Therefore, from now on 
we suppose q — I. By our notation we have 

i i—l 

Pr(^G, > r )= J>r(Y = j), i = l^...,l. (64) 

3=1 3=0 

We now bound I — E[Y] from above. 

i i i 

I - E[Y] = Z - E Pr(Y >i) = - Pr(E > i)) = E Pr <T < *) 

2=1 2 = 1 2 = 1 

/ 2 /.\ I i ^ i 

<=> E^(E G ,>-)sEPr(E G ,>E^) 

1=1 j = l j=l j = l J = l W 

< i + i>(iS>-^)i> E h 

i=l 3=1 ^ 3=i+l ^ J 



where (a) follows from (64); and (b) follows from Yl\=i ^-lv% < r - Now, by using Chebyshev's 



inequality 

Pr d E^ - > E -) ^ min ( X ' : ! , Vi G N, 1 < i < I - 1, (65) 

3=1 M 3=1+1 ^ VZ^j=i+l Pj J 

where var(Y^ =1 Gj) is the variance of the random variable YTj=i G r Therefore, 



i-i 



var (ELi Gj 



Z-£[E]<1 + Eniin(l, , 

i=l vZ_/3'=i+l pj ) 
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Due to independence of G;'s, we have 

Va r(^G,) = ^var(G J ) = ^^^<^^, < = 1, 2 I. (66) 

Therefore, 



v ■ *■ — 4 v - 



E*=i t 

/ - E'fy] < 1 + X)min(l, /" P ' ). (67) 

i=l \/-^j=i+lpj) 

But since pi > p 2 > . . . > we have Yj)=\ jj < ^ and (X5=;+i ^) 2 > (^f ) 2 - Therefore, 

i-i J_ z-i 

/ - E[Y) < 1 + ^min(l, ^) = 1 + ^min(l, jr^)- (68) 

i=i —pf~ i=i \ > 

Hence, by ( [68] ) and applying Lemma [6| which is stated and proved as following, the proof of 
Lemma [3] will be complete. 

Lemma 6: Assume I 6 N, and / > 1. Then, 



i-i 



1 + E min ( 1 '77^2)< Y + i- (69) 

i=l ^ ' 

Proof: For / < 18 the statement of the Lemma can be verified numerically. Therefore, 
suppose that / > 18. Let f(i) = , for i E N,l < i < I — 1; and consider the following three 
observations regarding the function /(.): 

1) f(i) increases as i increases for i e N, 1 < i < I — 1. 

2) /(I) = < I- 

3) /(/ - 1) = i=l > 1 

Therefore, 3m 6 N, 1< m < 1 - 1 such that 

m m + 1 

< 1 < vr—r — nrvr?- (70) 



(/-m) 2 - (/-(m + 1)) 2 



Note that m cannot be equal to I — 1, because ^_' ; _\^ 2 > 1- Rewriting the inequalities in (70) 
we get 



/- \!l+\-\<rn<l- Jl+\ + l- (VI) 



Now, we rewrite 1 + E!=i mm (l 



(l-i) 2 

l + ^min(l,— ^) = l + ^min(l,— -^) + ^ min(l, — -^-). (72) 



2-1 m 2-1 



(l-i)t J Z-^ v 'fZ-il 2 ^ (7-i) 2 

i=l v ; i=l v ; i=m+l v ; 
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But from (70) and the fact that f(i) 

l—l . m 



increases by increase of i, we have 



1-1 



1 + min(l, 



i=l 



{I -if 
l-l 



+ (/ — ! — m) = I — m + 



j=l—rn 



I- J 



l-l 



<l ' m+ ^ i(i-l) 

(a) m m 

< / — m + 



/ — m 



+ £ 



/ — m — 1 / 



-m+l 
2 



]=l-m 

I — m + 



j - j _ j - j 

J 1 J ' 



/ — m + 



/ — m 



i-i 

-- y - 
i ^ j 



j=l-m 



m(m + 1) 



{l-m - \){l-m + l- 1) 



(6) 
< 



3n 
4 2' 



51-9, I 



i i ii 



3' 

4 2' 



2\/l + - + 



-iJi + \ + ^i-eJi + \ + 3 



(73) 



where (a) follows from the Cauchy-Schwarz inequality; and (b) follows from f7T[ ). For I > 18 

6 



the term 



-V'+|+"'-6y'+|+3 . 



(V'+WX'+V'+i-i) 
true for all / > 1, / G N. 



in (73) is less than zero. Therefore, the statement of Lemma 



is 



Appendix E 
Proof of Corollary [T] 

Let LT* denote the partition (assignment) chosen by the optimal greedy static scheduling policy 
^-static- Therefore, we have ||-R(^. stat i c )||i = C T 3. Furthermore, consider an assignment, denoted 
by Udeu which maximizes the objective function in (JiJ). Let 77gfstatic denote the greedy static 
scheduling policy which corresponds to Tldet- Further, let ||i2det (^static) ||i designate the maximum 
number of objects that can be packed in the RP in ([5]) when a static scheduling policy ?7 st atic 

is implemented. Therefore, ||#det(^gfstatic)lli = c ^t, since 1 1 -Rdet (^gfstatic ) 1 1 1 is me value of 
objective function in (|5j) when the assignment is dictated by rj^^. The right part of the inequality 
in Corollary [T] in (jT]) is trivial since C T 3 is the optimal T 3 achievable under any scheduling policy. 
So we only need to prove the left part of the inequality in ([7]). Using a similar argument as the 
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one in part B of Section [TV] and by applying Cauchy-Schwarz inequality, we get 



ll^g-LtJIli > WR^irgLJWi - 2^/^(11^(^)111 + ^). (74) 



Now consider the function g(.) defined as follows: g(x) = x — 2y (N(x + ^)), xGl. 
So, g{x) is strictly increasing for x > ™. On the other hand, we know that 

C det = H^det^ltiJIll > I l^g* static) 111 > 1 1 %*- sta ti C ) 1 1 1 - N = C j3 - N, (75) 

where the right inequality follows from Theorem [T] By g(x) being an increasing function of x 
and f75j) we get 



I^hC/^u,,,)!!!. - 2 X /N(\ \R det (vl-LJ\\i + j)>C j3 -N- 2JN(C J3 - (76) 



Hence, by (74) and (76) we get 



ll%*Ltic)lii >C J3 -N- 2^N(C J3 - (77) 

Appendix F 
Proof of Theorem [4] 

By the same argument as the one in proof of Lemma [Tj C W _ T 3 can be achieved by a static 
scheduling policy in which the order of a client Rxj is according to the value of uijPij when 
assigned to APj. Therefore, by LLN, in order to achieve C W _ T 3, it is sufficient to find the 
assignment which has the highest expected weighted delivery for one interval. We start by 



proving the right part of the inequality in (48) 



First, we prove that for a given assignment II = [Xi,X 2 , . . . ,X N ] the optimal ordering of the 
packet of clients assigned to APj is according to the order of tOjPij, j E X;. To do so, it is 
sufficient to prove that for any given order of the clients if we swap the order of two adjacent 
clients in such a way that the client with higher UjPj is prioritized higher, then the expected 
weighted delivery will be no less than before swapping. 

Lemma 7: Let r, q E N, and ui, U2, ■ ■ ■ , co q E IKL Also, for some d E {1, 2, . . . , q — 1}, let 
uj[ = uji, for 1 < i < d and d+ 1 < i < q; and uj' d = u d+1 , u' d+1 = u d . Further, let Gi, G 2 , ■ ■ ■ , G q 
be independent geometric random variables with parameters pi,p 2 , ■ ■ ■ ,p q , respectively. Suppose 
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that uJdVd < ^d+iPd+i- In addition, let G'^G^, ■ ■ ■ , G' q be independent geometric random vari- 
ables, independent of GVs, with parameters pi,p2, • • Pd-i,Pd+iiPd,Pd+2, ■ ■ ■ ,p q , respectively. 
Then, 

q i q i 

1=1 j = l j=l jr' = l 

Proof: Let A = £Li co t Pr(£; =1 ^ < r), and = ELi w* Pr(Ej=i ^ < r). Then, 

B ~ A = ]T < < r) + J d Pr(£ < r) + Pr(]T G$ < r) 

i=l j=i j=i j=l 

<? j cZ — 1 i d 

+ E ^ Pr £ G ^ < ^) - E ^ Pr (E ^ < ^) - ^ *mE ^ ^ t ) 

i=d+l j=l i=l i=l i=l 

- u d+1 Pr(E < r) - Pr (E G ^ ^ r ) 

j=l i=d+l j=l 

d d+1 d 

= J d Pr(E G' 3 <r)+ J d+1 Pr(E Gj < r) - ^ Pr(E G 3 < r) 
j=i j=i j=i 

d+l 

j'=i 

= E Pr E <*i = *) K Pr ^ < ^ - + p r(G' d + G' d+1 < r - f )] 
t=i j=i 

- E Pr (E G i = *) Pr ( Grf - r - *) + Fv ( Gd + - r - *)] 
*=i j=i 

= E Pr E ^ = *) K Pr ^ < ^ - *) + p K^ + G' d+1 < r - 
t=i j=i 

- w d Pr(G d < r - t) - u d+1 Pr(G d + G d+1 < r - *)]■ 
Therefore, it is sufficient to show that for all t e N, 

^ Pr(G' d < t) + Pr(G' d + G' d+1 <t)-u d Pr(G d < t) - u d+1 Pr(G d + G d+ i < t) > 0. 

Note that 

. = w d+ i, and c^ +1 = w d . 



Pr(G" d < t) = 1 - (1 - and Pr(G, < t) = 1 - (1 - p d )*. 

+i)*-Pd+ 



Pr(G d + G d+ i < t) = Pr(G^ + G' d+l < t) = 1 - Pd(1 ' Pd+l)t ' Pd+l(1 ' Pd)t 
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Therefore, 

J d Pr(G' d <t) + u' d+1 Pr(G" d + G' d+1 <t)-u d Pr(G d < t) - co d+1 Pr(G d + G d+1 < t) (79) 

= (u d+lPd+ i - u dPd )( {1 - Pd+l)t - (1 - Pd)t ) > 0, teN, (80) 

Pd - Pd+i 

where the inequality follows from the assumption that u d+ ip d+ i — cu d p d > 0. ■ 



A. Proof of C W _ T 3 < C w . det + Nu 

r, 



We follow the same line of proof as in Section IV Since by the same argument made in the 
proof of Lemma [T| C W _ T 3 can be achieved using a static scheduling policy which uses the optimal 
ordering of clients according to ujjpfs, it is sufficient to show that for any static sheduling policy 
^wg-static which uses its corresponding optimal ordering we have 

^-T 3 (r] W g-static) < Ciu-det + Nu max . (81) 

Suppose an arbitrary static scheduling policy r/ W g-static with the corresponding partition n w g-static = 



[Xi,X 2 , . . . ,Xjv], which uses the optimal ordering is implemented. By (45) we know that 

M 

W-T 3 (r/wg-static) = WjRj (f/wg-static) • (82) 

On the other hand, by ([T]) we have 

Rj (Static) = liminf ^^(Mwg-stati^ J £ {1,2,..., M}. (83) 

r— s>oo T 

Let Yi denote the random variable indicating the total number of successful deliveries by A Pi 
during one interval, when 77 wg - s tatic is implemented, i = 1, 2, . . . , N. More precisely, 

y i = Yl N ^ ^g-static), i = 1, 2, . . . , N. (84) 

Define g« = |Xj|. Denote the enumeration of clients assigned to APj by {Xj(l),Xj(2), . . . ,Xj(gj)}, 
where the enumeration is according to the optimal ordering for the weighted case. Since a static 
scheduling policy is implemented and channels are i.i.d over time, by LLN we have 

ifrfa-)(w«rtc) = limi * f gkiW^^ = Pl{Yl > 3 ), i<j< ft , i<z<iv. 

r— >oo r 
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Therefore, 

N N Qi 

w-T 3 (r/ wg . static ) = Y ^^'(^wg-static) = Y Y ^(i) Fr ( Yi - 3) 
i=i jeXi i=i j'=i 

AT 9i j 

1=1 j=i fe=i 

Let Gij be a geometric random variable with parameter p^, z = 1, 2, . . . , N, j = 1, 2, . . . , M. 
Then, it is easy to see that 

k 

Yi = max k s.t. ^G iXi0) <r, ie{l,2,...,N},k< qi , (86) 

i=i 

since 77 wg . st atic persistently sends a packet until it is delivered, or the interval is over. The following 
lemma relates U and u/s to Yi. 

Lemma 8: Let 1 < u%, u 2 , ■ ■ ■ , u) q < co max for some co max E E. Also, let r G N and 
Gi, G 2 , ■ . ■ , G q be independent geometric random variables with parameters p 1 ,p 2 , . . . ,p q re- 
spectively, such that uipi > u 2 p 2 > • • • > u q p q > 0. Also define 

i 

I = max/ s.t. ^^l/j5j<r, (87) 



i=i 



and 



Then, we have 



F = maxi s.t. < r, i £ {1,2, ... , q}. (88) 

3=1 



q i I 

Y(J2 u i) Pr ( Y = t ) < Y UJ i + Uma *- (89) 

1=1 j = l jr' = l 

Hence, by Lemma [8] we have 

N qi 3 N h 

YYCY ^w) p*{Yi = j)<YY w ^o) + (90) 

i=l j=l k=l i=l j=l 



Therefore, by (85) and (90) we have 

N h 

w-t 3 (i] wg ., tatic ) <YY ux iti) + Nuj ™*- ( 91 ) 



1=1 j=i 



^ ut ' ^1=1X^^=1^0') * s me vaiue °f the objective function in (47) for a feasible solution. 
Therefore, 

N h 

Y Y - c ™-tef ( 92 ) 
i=i j=i 
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Putting ( |9T| ) and ( [92] ) together we get 

W-T (^wg-static) < Ciu-det + Nu max . (93) 

Hence the proof of the right inequality in Theorem [4] is complete. We now prove Lemma [8} 

Proof: Suppose that / > (for / = the proof is straightforward). We have 

i 1 J+i , 



i=i Pi «~T ^ 

Without loss of generality we can omit pj's that are equal to zero and assume < pi,p 2 , . . . ,p q < 
1. Furthermore, according to the same argument as in proof of Theorem 1, it is sufficient to 
prove the lemma for the case of q = r. Let Xi = Yl\=i ^i, where G{ = Geom(pi). We have 

t i l—l i 

J>) Pr(K = z) = J>) Pr(K = z) 

2 = 1 j=l j=l j=l 

T T 1 

+ ^Pr(X l = t)(X;(E w i) Pr ( y = i l^ = *)) 

t=l i=l j=l 

However, since u max > uipi > coi+m+i > ■■■ > u) T p T > 0, J2i=i(J2)=i u i) Pr (^ = is less 
than the case where uipi = ui + ipi + \ = . . . = co T p T . With a similar argument as in the proof of 
Theorem 1 we get 

t i l—l i r I 

EE ^) Pr ( r = **) ^ EE ^ Pr ( y = o + EE ^ + - £ )) Pr ( x < = *) 

i=l j=l j=l j=l t=l j=l 

l-l i I 

= EE^) Pr ( y = o + E w i + - Pr ^ > r )) 

i=i j=i j=i 

oo oo (— 1 i I 

-^E tPr ( x ' = t )- E ^pr(^ = t)] = EE^) pr ( r = z ) + E^+^ r ) 

t=l t=T+l 1 = 1 j = l j = l 

I l—l I oo 

- e + ujipit) J2 Pr ( y = *) - "ipi E ~ + uipi E * Pr ( x ' = *) 

j=l 1=0 i=l ^ l t=r+l 

l-l i I I 1 I 

= EE w j ~ E u i ~ uipiT ^ Pr ( r = ^ + E w j + uipi ( T ~ E ~)) 

1=0 j=l j=l j=l 8=1 ^* 

OO 

+ W(p, ^ tPr(X, =t) 

t=T+l 

^ l — l i I I oo 

< EE w j ~ E u o ~ uipiT "> Pr ( y = i ') + E ^ + + w ^ E 1 Pr ( x * = *) 

8=0 j = l j=l j = l t=T+l 
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l-l I I l-l I 

-£(- £ wi-^)Pr(y=i) + (Ew i +w I+1 )+x;p^=o(w J p l r+ £ ^; 

i=0 j=i+l j=l i=o i=«+i J 



i ; — i ; i ,j i z 

E w j + + £( £ "f^ ~ E w j) Pr ( y = *) - E u i + - E w * + 

j=l i=0 i=i+l ^ j=i+l j=l j=l 



where (a) follows from r — X)i=i ^ < an ^ — ^i+iPi+i', (b) follows from (62); and (c) 
follows from the fact that Vj e {i + 1, . . . , 1} ^ < w 7 . ■ 

B. Proof of C w . det - 2u max JN(C w . det + f ) < C TO _ T 3 

The proof of the lower bound is similar to the proof of lower bound in Theorem [TJ Consider 



the assignment proposed by the solution to the problem in (47), where the clients which have 



not been assigned to any AP for transmission are assigned to AP's arbitrarily. Let II^L 



g-static 



[Xf et ,X2 et , . . . ,X^ 1 ] denote the resulting partition, and also let ^g-statk denote the correspond- 
ing static scheduling policy which orders clients based on their channel success probabilities. 



Therefore, w-J 3 (r)^_ sMic ) < C W _ T 3. So, it is sufficient to prove that C w . det -2^ N(C w . det + f ) < 
w '~^ (^wg-static)- Let Wf &t denote the random variable indicating the total weight of successful 
deliveries by APj during one interval, when ^g-static * s implemented, i = 1,2, . . . , N. More 
precisely, 

Wi = "MO-* ^-static), i = 1, 2, . . . , N. (95) 

j£Zf< 

Then, by LLN we have 

TV 

^-T 3 « t . statlc ) = ^i?[^ t ]. 

i=l 

Therefore, it is sufficient to prove that 

I N 

C w . det - 2JN(C w . det +-)< E W?\ 

i=i 

Define qi = |Zf et |, and denote the enumeration of clients assigned to APj by {Zf et (l),Zf et (2), . . . 
,Zf et (gj)}, where the enumeration is according to the channel success probabilities of differ- 
ent clients in Zj. Further, let Gij be a geometric random variable with parameter p^, i = 
1, 2, . . . , N, j = 1, 2, . . . , M. It is easy to see that 

k k 

Wf* = max^u) j s.t. ^ G^-* < r, i e {1, 2, . . . , N}, k < Qi , (96) 

i=i j=i 
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since 77!^^ persistently sends a packet until it is delivered, or the interval is over. Also define 



Zf et = max I s.t. \J < r, Z < g$. 

i=i 

Therefore, Zf et is the maximum number of objects that fit into a bin of capacity r when the 
channels are relaxed and clients in lf et are assigned to AP^. By Lemma [9j which is stated and 
proved later in this appendix, we have 

N I N 



i=\ j=i 



E E ^ - aw E A det + \ < E E i w " st ] 



i=l 



(97) 



1=1 



Note that 5X1 E*= 1 — Cw-det- According to Cauchy-Schwarz ineqality we have 



i=l 



UJ r , 



A 



/V 



iV(^Zf + 



i=l 



< Wr. 



TV If 



, ^(E E Ux ^) + = N ( c *>** + t% 



=1 i=i 



(98) 



Putting (97) and (98) together we get 



N 

C w -det — 2y iV(C to .det + — ) < C„,.T 3 - 



(99) 



Hence, the left inequality of Theorem [4] is proved and the proof of Theorem [4] is complete. We 
now state and prove Lemma |9j 

Lemma 9: Let 1 < ui\, ■ ■ ■ , oj q < u> max for some u) max G R. Also, let r G N and 
Gi, G 2 , • ■ ■ > G q be independent geometric random variables with parameters p 1 ,p 2 , . . . ,p q re- 
spectively, such that 1 > pi > p 2 > . . . > p q > 0. Also define 



Z = maxZ s.t. 1/Pi — r ) 



and 



Then, we have 



K = maxi s.t. E^ Gj < r, i G {1, 2, . . . , q}. 

3=1 



/ I 91 



(100) 



(101) 



(102) 



i=i j=i 
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Proof: With the same argument as in Theorem 1, it is sufficient to assume q = I. The proof 

is very similar to the proof of lower bound in Theorem 1 : 

i i i ill 

5> - E(E w i) Pr ( r = o = 5> - 5>(E Pr ( y = j)) 

i=l i=l j'=l i=l «=1 j'=i 

= ^o; l Pr(G 1 + G 2 + ...G J >r) 



i=i 

<E^ Pr (iE(^--)i> E -)+^<E^ min (^ t^i M 



i-i 4 n 

< ^ mm(l, — - x — ) + w, < wj + 2^ Wi mm(l, ) 

i=i \2^j=i+i Pj J i=i \ ' 

i-i 



max \/ 1 i . ■ 



i=l 



where (a) follows from Chebyshev's inequality; and (b) follows from Lemma [6j 
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